Blume-Emery-Griffiths Model in a Random Crystal FieldB 
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We study the Blume-Emery-Griffiths model in a random crystal field in two and three dimensions, 
through a real-space renormalization-group approach and a mean-field approximation, respectively. 
According to the two-dimensional renormalization-group calculation, non-symmetry-breaking first- 
order phase transitions are eliminated and symmetry-breaking discontinuous transitions are replaced 
by continuous ones, when disorder is introduced. On the other hand, the mean-field calculation 
predicts that first-order transitions are not eliminated by disorder, although some changes are in- 
troduced in the phase diagrams. We make some comments on the consequences of a degeneracy 
parameter, which may be relevant in martensitic transitions. 
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I. INTRODUCTION 

The Blume-Emery-Griffiths (BEG) model is a spin- 
f Ising model, originally proposed to study 3 He- 4 He 
mixturestl. Later, it was used in the description of a va- 
riety x>f different physical phenomena: multicomponent 
fluidsa, microemulsionscL semiconductor alloysu, elec- 
tronic conduction modelsD, etc. Its hamiltonian reads: 

n = - J E s * s s - K J2 s * s j + E A ^ 



<»,i> 



<*,]> 



where the first two sums are over all nearest-neighbor 
pairs on a lattice, the last one is over all sites and Si = 
±1, 0. J is the exchange parameter, K is the biquadratic 
interaction and Aj is a site-dependent crystal field (A, = 
A for the pure model) . The phase diagram of the model 
presents first-order and continuous phase transitjans and, 
for K < 0, a rich variety of multicritical pointsail 

Nevertheless, some systems were found to be better 
described by a disordered BEG model, as, for instance, 
conventional shape memory alloysB and fluid mixtures on 
disorder materials (like aerogel, for example p^3. From 
the theoretical point of view, on the other hand, it has 
been argued that randomness may ..have drastic conse- 
quences on multicritical behaviortilE3. In two dimen- 
sions, for instance, any infinitesimal amount of disor- 
der supresses non-symmetry-breaking first-order phase 
transitions and replaces symmetry-breaking first-order 
phase transitions by continuous ones. The effect of disor- 
der on three-dimensional systems is different: first-order 
phase transitions only disappear at a finite amount of 
randortmcp i3. This behavior has been observed in some 
models£3~E£ . 

In order to study the effects of disorder on its phase- 
diagram, we study the BEG model in a random crystal 
field (henceforth called RBEG model) given by the prob- 
ability distribution: 



P(A,-) = r 6(A Z + A) + (1 — r) 5(A< - A) 



(2) 



ft is worthy stressing that the exact form of the disor- 
der is not relevant to the overall qualitative consequences 
on the phase diagram. If randomness is chosen to be 
in the interactions J or K, the qualitative effects will be 
the same (in what concerns first-order phase transitions). 
This is due to the fact that, even if the initial disorder 
is on the bonds (interactions J or K), a scale transfor- 
mation will propagate this disorder to the crystal field 
term, which will act just like field randomness on the 
coexistence boundary. Moreover, the exact form of the 
probability distribution is not relevant, either; we have 
performed calculations with other distributions and they 
lead to the same qualitative picture as the one found in 
this work. 

Finally, we would like to mention that, to the best 
of our knowledge, the BEG model in a random crystal 
field has not been studied so far. Previous studjpSjXaa.- 
centrated on the random Blume-Capel modelBEJO'ESEZl, 
which has a simpler phase-diagram than the BEG 
model's. 

The remainder of this paper is organized as follows. 
In section II we outline the mean-field approximation 
we use and discuss the results, in section III we present 
the real-space renormalization-group (RSRG) calculation 
(expected to hold for two-dimensional systems), and in 
the last section we summarize our main conclusions and 
comment on the influence of a degeneracy parameter p 
on the critical behavior. 



II. MEAN-FIELD CALCULATION 

We chose an ordinary mean-field approximation to 
study the three-dimensional system. The procedure is 
rather usual and we refer the reader to Ref. [ll| where a 
detailed discussion of the method is done. 
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However, we would like to stress that the mean- field 
approximation we use is equivalent to a model where the 
interaction is of infinite-range, i.e., each spin interacts 
with every other spin in the system. This will have ex- 
plicit consequences on the phase diagram and we will 
return to this point later. 

Most of the information about the phase diagram is 
numerically calculated but some analytical results can 
be obtained. Among them, we can find the ground state 
for any values of J, K, A and r. 

It is possible to show that the ground state magnetiza- 
tion, too, for d = A/zJ > 0, where z is the coordination 
number of the lattice, is given by [results for A < can 
be inferred from the mapping (r, A) <-> (1 — r, — A)]: 



to 



l-(l-r) 



d-(k+l) 



1 



(3) 



where k = K/J and 6[x] is the step function, such that 
9[x] = or 1 for x < or x > respectively. The 
ferromagnetic phase 0\ (see figures in this subsection), 
with mo = 1, is stable for d < d c = (k + 1) {^-), while 
for d > d c the ground state is such that too = r (denoted 
02 in our figures). Note that, except for r = 0, the 
ground state is always ordered; this is a consequence of 
the simple mean-field approximation we used (we will 
return to this point below). 

One can obtain the continuous transition line exactly, 
by expanding <& m in in powers of the magnetization to 
and taking the coefficient of to 2 equal to zero: 



(4) 



-kpdjt c 



2 + e~ k e 



-d/t c 



where t c = ksTc/zJ. More specifically, note that, for 
d > 1, the value of the critical temperature is t c = r. 
So, for any value of r ^ 0, the critical line between the 
paramagnet and the O2 phases extends to d = 00 (see 
figures in this subsection). This is not the expected be- 
havior for a cubic lattice, for the following reason. The 
RBEG model for d = 00 is equivalent to the site-diluted 
spin-1/2 Ising model, since for A = 00, a +A crystal 
field acting on a given site forces that site to be in the 
S = state (absent), while a — A field forces the site to 
be either in the state S = 1 or in the state S = — 1 (both 
represent a present site). Thus, only for high enough r 
an infinite cluster of S = ±1 states will form and will be 
able to sustain order. Exactly at r = r c , there is such an 
infinite cluster but its critical temperature is zero. There- 
fore, the critical parameter d c should only reach infinity 
for r > r c . However, the simple mean- field analysis we 
made leads to r c = 0, since it is equivalent ta .a- model 
with infinite-range interactions. In some casesaEZl, more 
elaborated mean-field- like procedures were applied to the 
Blumc-Capel model in a random crystal field. Briefly, 
the consequence of these approaches is that the transi- 
tion line between O2 and D phases does not extend to 
d = 00 for all values of r. All other results are similar 



to the ones obtained with our simple mean-field approx- 
imation. We note in advance that the approach we used 
for the two-dimenisonal model leads to a finite value of 
r c , as expected on physical grounds. 

We have already pointed out that t c (A = 00) does 
not depend on K\ this comes from the mapping be- 
tween the RDBEG model and the site-diluted spin-1/2 
Ising model. The S — states (absent sites) play no 
role in the dynamics of the model and the present sites 
can only be in the states S = 1 or S = — 1; thus, the 
biquadractic interaction, K, is irrelevant in this limit. 
If, for instance, the probability distribution 'P(Aj) = 
r 6(Ai) + (1 - r) 8(A l - A) is used, the A = 00 limit will 
be equivalent to the site-diluted spin-1 Ising model; then, 
t c (A = 00) will depend on K. Note that the discussion 
in this paragraph applies to the two-dimensional case as 
well. 

We now turn to the discussion of the ksT/zJ x A/zJ 
phase diagrams. In Figs. H, |2[ || and f| we depict sections 
of constant K/J — 5, for many values of r. The phase 
diagram for r — (pure BEG model) is qualitatively the 
same as for r — 0.1 (Fig. [I]), except that the O2 phase is 
not present. 
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FIG. 1. Mean-field phase diagram for K/J = 5 and r — 0.1. 
Filled circles and the open square stand for critical end points 
and a critical point, respectively. Continuous (dashed) lines 
represent continuous (first-order) transitions. The phases are: 
ordered with m = 1 (Oi), ordered with m — r (O2), disor- 
dered with q > 1/2 (-Di), and disordered with q < 1/2 (D2). 



Note that the size of the ordered phases increases with 
r. This is expected, since r is the fraction of sites which 
feel a —A crystal field (we have already commented on 
the "tail" which separates the O2 and D2 phases, given 
by t c = r). Another important feature is the presence 
of a first-order line between two disordered phases, for 
r = 0.1 and r = 0.3. In both of these phases to = but 
q > 1/2 for Di, while q < 1/2 for D 2 . We would like to 
call attention for the phase diagram for r = 0.3 (Fig. |^); 
this type of diagram is not present in the Blume-Capel 
model. 
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FIG. 2. Mean-field phase diagram for K/J = 5 and r = 0.3. 
Same conventions as in Fig. jjj; T stands for tricritical points. 
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FIG. 3. Mean-field phase diagram for K/J = 5 and r = 0.5. 
Same conventions as in Fig. |2| 



those of Refs. [lO] and and we shall not depict all of 
them here. The only exception is the diagram for r = 0.1 
(Fig. ||), which is not present for higher values of K/J. 
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FIG. 5. Mean-field phase diagram for K/J = and r = 0.1. 
Same conventions as in Fig. |2| 



On general grounds, one should note that the mean- 
field approximation we employed suggests that the ran- 
dom crystal field does not destroy the first order transi- 
tions between disordered phases and between an ordered 
and a disordered phase. Even first order lines between 
ordered phases (like the one in Fig. [sj) survive the intro- 
duction of randomness. 



III. TWO-DIMENSIONAL 
RENORMALIZATION-GROUP 
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FIG. 4. Mean-field phase diagram for K/J — 5 and r = 0.7. 
Same conventions as in Fig. |j| 

The corresponding phase diagrams for K/ J — 3 show 
only three types of critical behavior: for r near zero, they 
are similar to the phase diagram for K/J = 5 and r = 0.3 
(Fig. ||); for intermediate values of r, the behavior is the 
same as for Kj J = 5 and r = 0.5 (Fig. ||); and for r close 
to one, the equivalence is with the diagrams for K/ J = 5 
and r = Q.7 (Fig. %. 

The Blume-Capel model (K/ J = 0) has .already .been 
studied within mean-field approximationaJ't3E£fE3, al- 
though for different probability distributions; the results 
we find in this case are in qualitative agreement with 



It is well known that mean-ficld-like approximations 
are not suitable to describe low-dimensional systems. We 
have then to resort to a different technique, in order to 
study the RBEG model in two dimensions. RSRG proce- 
dures, on the other hand, have been successfully applied 
to two-dimensional systems. Note, however, that RSRG 
approximations, in general, do not lead to results as pre- 
cise as those obtained with Monte Carlo simulations, phe- 
nomenological renormalization or conformal invariance 
techniques. Nevertheless, they allow for a correct de- 
scription of universality classes, order of the transitions, 
crossover phenomena, etc. 

The procedure is the same as the one thouroughly dis- 
cussed in Ref. [l3| There is just one technical point we 
would like to stress. Although we start with a uniform 
distribution for J and K, the renormalization procedure 
will introduce randomness in all renormalized quantities 
(J', K 1 and A'). A possible approach is to follow the suc- 
cessive renormalized distributions of these parameters in 
order to study the phase diagram. We adopted an alter- 
native way, which forces the renormalized distributions to 
be the same as the initial ones, but with renormalized pa- 
rameters, namely, P' ap (J) = S(J - J'), V (K) = S(K - 
K 1 ) and V' ap (Ai) = r' S(Ai + A') + (1 - r') <$(A 2 - A'). 
The values of J' and K' are obtained imposing that the 
first moment of the actual distributions for J and K and 
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of V (J) and V' ap (K) are equal, respectively. The val- 
ues r and A' are calculated imposing that the two lowest 
moments of V' ap (A) match those of the real distribution. 
This procedure has to be used with some care: in certain 
systems, where the random-field mechanism is important 
and the initial randomness is on the interaction (J, for in- 
stance) , forcing the field back into a uniform distribution 
leads to incorrect results. In Ref. |l^, for instance, the 
crystal field probability distribution is maintained uni- 
form throughout the renormalization procedure. Con- 
sequently, the critical behavior of the random model is 
characteristic of a high-dimensional system: the critical 
temperature of the tricritical point diminishes as random- 
ness is increased but only reaches the zero temperature 
axis at a finite value of the disorder. As discussed in Ref. 
|l2| , the mechanism responsible for the lack of first-order 
phase transitions in two-dimensional random systems is 
the disorder in the crystal field, which is not taken into 
account by approximations such as the one used in Ref. 
|l9| . In the model we study in this paper, however, the 
disorder in the field is not approximated away by our 
RSRG procedure. 

Our results for K/ J = 2 are presented in Fig. ^, where 
we depict the kT/zJ x A/zJ phase diagram for r = 
(pure BEG model), r = 0.2, r = 0.45, and r = 0.5. 

Let us first comment on the pure BEG model (curve 
(a) of Fig. |^). As for K/J = 5 in three dimensions, 
there are two types of disordered phases: both have 
to = but q > 1/2 for phase D\ and q < 1/2 for 
phase D 2 . The continuous line between phases O and 
D\ belongs to the universality class of the Ising model: 
this line is attracted to the Ising fixed point, C* = 
(J = 0.4407, K = -0.07308, A = -co). The dashed line 
between phases O and D2 is attracted to the fixed 
point F\ = ( J = 00, K = 00, A = 2( J + K)), which rep- 
resents a first-order transition in both to and q, i.e., 
the largest eigenvalue of the even and the odd sec- 
tors of the linearized RGT matrix are equal to b d (see 
Ref. On the other hand, the dashed line between 

phases D\ and D2 is attracted to the fixed point F2 = 
( J = 0, K = 00, A = 2K + In 2); in this fixed point only 
the largest eigenvalue of the even sector of the linearized 
RGT matrix is equal to b d ; this is a sign of a discontinuity 
in q (but not in to) when the line is crossedE3. 
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FIG. 6. Renormaliza- 

tion-group phase diagram for K/J = 2 and (a) r = 0, (b) 
r = 0.2, (c) r — 0.45, and (d) r = 0.5. Filled circles stand for 
critical end points, O for the ordered phase and D\ and D2 
for the two disordered phases (see text). Continuous (dashed) 
lines represent continuous (first-order) transitions. The tran- 
sition lines extend to A — » 00 only for r > 0.5. Note that the 
critical end point (filled circle) is present only for r — 0. 

In curves (6), (c), and (d) of Fig. ^ we depict the 
kT/zJ x A/zJ phase diagram for r ^ 0. We note that 
the first-order line is either replaced by a line of contin- 
uous transitions (between O and D2 phases) or is elimi- 
nated (between D\ and D2 phases), for any infinitesimal 
amount of randomness. In fact, the first-order fixed point 
attractors, F\ and F2, are unstable along the r direction. 
There is still a line separating the two disordered phases 
(not depicted in Fig. ||), D\ and D 2 , for r^0, but this 
line is attracted to the (r = 1/2, J = 0, K = 0, A = 00) 
fixed point. This point represents a model with indepen- 
dent spins, in which no phase transition can take place. 
We note that our results are in accordance with general 
arguments on the effects of randomness on multicritical 
phase diagramsEJ, although, to the best of our knowl- 
edge, some features of these arguments have never been 
tested so far. 

On the other hand, the whole line of continuous tran- 
sitions for r / belongs to the pure Ising model uni- 
versality class, i.e., C* is a stable fixed-point along the 
r direction. This is expected, since, for the hierarchical 
lattice used in this work, the specific heat critical expo- 
nent of the pure Ising model, a, is negative and disorder 
is irrelevant, according to the Harris criterioncil. For the 
corresponding model on a two-dimensional Bravais lat- 
tice, where a = 0, the Harris criterion is inconclusive. 
The accepted behavior, when disorder is present, is the 
following: critical exponents of the random model retain 
the same values as their pure conterparts but logarithmic 
corrections are introduced by randomnessE3. Experimen- 
tal results also indicate the same critical exponents for 
pure and random two-dimensional Ising modelo. 

We would like to call attention to the behavior of the 
critical point which separates the O and D2 phases at 
T = 0. For r < 0.5, the transition at zero tempera- 
ture takes place at a finite value of A/zJ. For r > 0.5, 
the critical line between the ordered and the disordered 
phases extends to A/zJ — 00 in the diagram. In fact, for 
A/zJ = 00 the RBEG model is equivalent to the site- 
dilute spin- 1/2 Ising model, as discussed above. Thus, 
only for high enough r an infinite cluster of S = ±1 states 
will form and will be able to sustain order. There is a crit- 
ical value of r, r c , such that the critical line between the 
ordered and disordered phases only reaches A/zJ = 00 
for r > r c . Our evaluation of r c is 1/2, while the accepted 
value for the site percolatiotucritical probability on the 
square lattice is r c = 0.5927EJ. This difference is due to 
the small cell we use in this work; nevertheless, the cor- 
rect qualitative picture is obtained, i.e., a finite value of 
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r c - 

Finally, we would like to stress that there are only two 
types of phase diagrams for the BEG model; for high 
values of K/J these diagrams have the same structure as 
for K/J — 2. For small values of K/J, the structure is 
the same as for the Blume-Capel (K, = 0) model. As this 
model has been studied elsewherd!3, we will not discuss 
it here. 



IV. SUMMARY 

We studied the BEG model in two and three dimen- 
sions within a RSRG framework and a mean-field ap- 
proximation, respectively. The disorder is on the crystal 
field term, which follows a probability distribution given 
by: V(Ai) = r <5(A; + A) + (1 — r) <5(A; - A). 

For the mean-field approximation (expected to repre- 
sent the qualitative behavior of three-dimensional sys- 
tems) , the presence of randomness increases the ordered 
phase and brings qualitative changes to the kT/zJ x 
A/zJ phase diagram. More specifically, first-order tran- 
sitions are present in the disordered model, but new mul- 
ticritical points emerge, depending on the value of r. 

In two dimensions, the RSRG approach we use shows 
that randomness has a drastic effect on critical behav- 
ior: it supresses non-symmetry-brcaking first-order tran- 
sitions and replaces symmetry-breaking discontinuous 
transitions by continuous ones. These results are in ac- 
cordance with general arguments concerning the effects 
of quenched impurities on multicritical behavior (but, to 
the best of our knowledge, the disappearance of the first- 
order line between disordered phases or between ordered 
phases has never been seen in an actual calculation). 
The line of continuous transition, present for the disor- 
dered (r ^ 0) model, belongs to the universality class of 
the two-dimensional pure Ising model; this results agrees 
with the Harris criterion, since the specific heat critical 
exponent, a, is negative for the hierarchical lattice used 
in this work. It has been conjectured that a new unsta- 
ble critical point, at finite .temperatures, might be present 
for the disordered systemH; we found no evidence of this 
point, for any value of K. 

We have also studied the so-called degenerate Blume- 
Emery-Griffiths (DBEG) model, introduced in the study 
of martensitic transitions^. In the DBEG model, the 
S = states are assumed to have a degeneracy p, which 
mimics the effects of vibrational degrees of freedom. It 
has been shown in Ref. || that the effect of increasing p 
is to shrink the ordered phase and to increase the re- 
gion where the transition is of first-order. Using the 
same probability distribution for the crystal field as in 
the RBEG model, we were able to show that the param- 
eter p may bring only quantitative changes to the phase 
diagrams, for any K/J, r, and p. This is expected, since 
the DBEG model is equivalent to the usual BEG model 
with all crystal fields displaced by ln(p)El In particular, 



any infinitesimal amount of randomness in two dimen- 
sions destroys the first order transitions, no matter the 
value of p. 

Finally, we would like to stress that our approxima- 
tion does not allow for a study of the BEG model with 
negativa.iC, where new and interesting critical behavior 
emergesQ. Work is now being made to discuss this model 
in the presence of a random crystal field. 
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